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Abstract 

The g-state Potts field theory describes the universality class associated to the spontaneous 
breaking of the permutation symmetry of q colors. In two dimensions it is defined up to 
q — 4 and exhibits duality and integrability away from critical temperature in absence of 
magnetic field. We show how, when a magnetic field is switched on, it provides the simplest 
model of confinement allowing for both mesons and baryons. Deconfined quarks (kinks) exist 
in a phase bounded by a first order transition on one side, and a second order transition 
on the other. The evolution of the mass spectrum with temperature and magnetic field is 
discussed. 



1 Introduction 



Confinement is that property of quantum field theory for which excitations which are genuine 
asymptotic particles in a region of coupling space become unobservable in isolation in another 
region, where they leave the place to new asymptotic particles (mesons, baryons, ...) of which 
they can be seen as "constituents" (quarks). The case of two-dimensional space-time provides 
the framework in which this, as other general properties of quantum field theory, can be studied 
in their simplest form. It is known that the gauge theory setting in which confinement is 
ordinarily discussed in four dimensions becomes somehow redundant in two dimensions. Indeed, 
due to the absence of transverse spatial dimensions, massless gauge fields do not carry particle 
degrees of freedom in d = 2, so that an alternative description of the theory exists which relies 
on physical excitations only. For example, two-dimensional quantum electrodynamics can be 
exactly mapped onto the theory of a self-interacting neutral boson which makes quite transparent 
the presence of quark confinement [1, 2]. In absence of electromagnetic interaction the quarks 
correspond to the solitons interpolating between the vacua of a periodic bosonic potential. Quark 
interaction destroys the degeneracy of the bosonic vacua and removes the topologically charged 
excitations from the spectrum of asymptotic states. What remains is a spectrum of mesons 
originating from confinement of soliton-antisoliton pairs 1 . 

This mechanism of confinement through breaking of degeneracy of discrete vacua is quite 
general in two dimensions and exhibits its most essential features in the case of a finite number 
of vacua originating from the spontaneous breaking of a discrete symmetry, with the kinks 
interpolating between these vacua playing the role of the quarks. Then it is not surprising that 
two-dimensional Ising field theory (i.e. the field theory describing the scaling limit of the two- 
dimensional Ising model) provides the simplest model of confinement (only two vacua). The 
associated mesonic spectrum was first studied in [3]. 

In theories, like Ising, with a one-component order parameter, the confined particles are 
made of an even number of quarks. Indeed, in this case the vacua are located along a line in 
order parameter space 2 , so that the kink sequences starting from and going back to the true 
vacuum (the only ones generating bound states via confinement) consist of a number 2j of kinks. 
The lightest bound states are mesons (kink-antikink composites) corresponding to j = 1. 

On the other hand, when the order parameter has more than one component we can find 
three vacua located on a plane in order parameter space. Now we can have a three-kink sequence 
making a loop through the three vacua. Such a sequence is confined into baryons, a finite number 
of which must be stable sufficiently close to the deconfining point. Indeed, if m is the mass of 
the kink, the lightest baryons will have mass ~ 3m and will not be able to decay into two mesons 
with mass ~ 2m each. 

1 We refer to the generic case in which the 0-angle is not fine-tuned to the specific value which partially preserves 
vacuum degeneracy. 

2 This applies also to two-dimensional quantum electrodynamics, with the notion of one-dimensional order 
parameter referred to the bosonic version. See [4] for a discussion of confinement in theories with more general 
bosonic potentials. 
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In this paper we consider the simplest model of confinement allowing for baryons, i.e. the 
field theory describing the scaling limit of the two-dimensional g-state Potts model. The latter 
generalizes the Ising model to the case in which each site on the lattice can take q different 
colors, and has an order parameter with q — 1 components [5, 6]. In absence of magnetic field, 
the ferromagnetic g-state Potts model undergoes an ordering transition which is continuous up 
to a number of colors q c , which in two dimensions equals 4 [7]. The scaling limit then can be 
taken up to this value of q and produces a field theory in which a magnetic field confines the 
kinks and allows for baryons at q = 3, 4. 

In absence of the confining field, the quarks (kinks) behave as free neutral fermions in the 
Ising case (q = 2), but become interacting for q = 3,4. For q = 4 the zero-field theory is 
equivalent to the sine-Gordon model at a specific value of the coupling (/3 2 = 2ir, see e.g. [8]), 
while for q = 3 the ultraviolet fixed point is non-trivial. In any case the confining interaction 
is non-local with respect to the quarks. The possibility of a quantitative study in weak field 
comes from integrability of the g-state Potts field theory in zero-field [9]. Hence, form factor 
perturbation theory [10] allows to express mass corrections in terms of the matrix elements of 
the magnetic operator computed in [11]. 

Here, however, our main interest will be in a qualitative characterization of the evolution 
of the mass spectrum for generic values of the temperature and of a magnetic field chosen to 
act on a single color. For q = 3, 4 this choice allows for an extended phase on the parameter 
plane in which the quarks are deconfined. Such a phase is bounded by the confining (first order) 
transition on one side, and by a spontaneous breaking (second order) transition on the other. 
Outside this region, the spectrum of asymptotic particles is made of mesons (everywhere) and 
baryons (at least sufficiently close to the deconfining transition). 

The paper is organized as follows. In the next section we discuss the g-state Potts model 
with our choice of magnetic field, starting with the lattice definition and then switching to the 
field theoretical description of the scaling limit. In section 3 we focus on the two-dimensional 
case and recall the exact scattering solution for the low-temperature phase in zero-field, before 
showing how this is related by duality to the scattering solution for the high-temperature phase. 
Section 4 is devoted to the weak field analysis, while the spectrum evolution as a function of 
temperature and magnetic field is discussed in section 5. Few final remarks are collected in 
section 6. 

2 Potts model with magnetic field 

In this section we discuss the lattice definition of the g-state Potts model and the field theoretical 
description of the scaling limit. Some remarks about d > 2, and about q > 4 in d = 2 are included 
although they are not used in the rest of the paper. 
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2.1 Lattice model 

The (/-state Potts model [5, 6] is a generalization of the Ising model in which each site variable 
s(x) at site x on the lattice can assume q different values (colors). In absence of magnetic field 
the interaction only distinguishes whether nearest neighbor sites have equal or different color, 
so that the Hamiltonian is invariant under the group S q of permutation of the colors. If we 
add a magnetic field H acting only on the sites with a specific color (say s = q), the reduced 
Hamiltonian can be written as 

and is invariant under the group S q -i of permutations of the first q — 1 colors (the first sum is 
over nearest neighbors). 

In the ferromagnetic case at H = 0, the q configurations in which all the sites have the same 
color minimize the energy and the system exhibits spontaneous magnetisation for sufficiently low 
values of the temperature T. Above a critical temperature T c the thermal fluctuations become 
dominant and the system is in a disordered phase. If we introduce the variables 

°~a(x) = 8 s (x),a ~ ~ , a = l,2,...,q (2) 

satisfying the condition 

Q 

£> Q (x)=0, (3) 

a=l 

the expectation values (a a ) differ from zero only in the low-temperature phase and can be used 
as order parameters. 

When the magnetic field is switched on with a positive value, the ground state at T = 
is unique (all sites have color q), and there can be no phase transition as the temperature is 
increased. 

Different is the situation for H < 0, q > 2. As H — > — oo the color q becomes forbidden and 
a zero-field (q — l)-Potts model is obtained. Then the critical points at H = and H = — oo are 
the endpoints of a phase transition line which in the T-H plane separates a low-temperature, 
spontaneously magnetized phase with q — 1 degenerate ground states from a high-temperature, 
disordered phase (Fig. 1). 

The nature of the transition at T = T C (H), H < 0, depends on q and on the dimensionality 
d. It is well known (see [6]) that in the zero- field g-state Potts model there exists a value q c (d) 
(not necessarily integer 3 ) such that the transition is continuous for q < q c and first order for 
q > q c . Accordingly, three cases can be distinguished for the transition going from C q to C q -\ 
in Fig. 1: 

3 One can make sense of the Potts model for non-integer values of q through the mapping onto the random 
cluster model [12]. Although we will have mainly in mind integer values, most of our discussion can be done 
treating q as a continuous parameter, and this will be understood in the following. 
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Figure 1: Phase diagram of the model (1) for q > 2. The values H = and H = — oo correspond 
to the zero-field q- and (q — l)-state Potts model, respectively. The ordered phase possesses q—l 
degenerate ground states and is separated from the disordered region by two phase transition 
lines (dashed and thick line) . The dashed line is a first order transition line along which q ground 
states are degenerate. The nature of the transition along the thick line depends on q and d. 



i) q — 1 > q c . The transition is first order with q phases coexisting along the transition line. 
Only at C q q + 1 phases coexist. 

ii) 2 < q < q c . The transition is continuous. 

iii) q — 1 < q c < q. The correlation length is finite at C q and infinite at C q ~\. The nature of the 
transition depends on d. 

The value f/ c (2) is exactly known to be 4 [7], while (/ c (3) lies in between 2 and 3. Hence, in 
d = 2 the transition induced by the field is continuous for q = 3, 4 and first order for q > 5; in 
d = 3 it is first order for q > 4. The cases q = 5 in d = 2 and q = 3 in d = 3 are of the type iii) 
above and will be discussed in a moment. 



2.2 Field theory description 

A continuous, field theoretical description at scales much larger than the lattice spacing is 
possible at and around those points of the phase diagram where the correlation length diverges. 
For q < q c , the transition point at H = {C q in Fig. 1) corresponds in the scaling limit to a 
fixed point of the renormalization group, i.e. to a conformal field theory with action A^ FT . The 
scaling limit of (1) around C q is described by the action 

A = Aqpt -r J d d x e(x) - h J d d x a q {x) , (4) 

where e(x) is the leading 5 g -invariant operator in A^) FT , and and a q {x) is the leading S q ^\- 
preserving magnetic operator. If X$ denotes the scaling dimension of an operator $(x) at the 
S^-invariant fixed point and m a mass scale, the couplings r and h behave dimensionally as 

r - m d ~ x - q) , h ~ m d - x ° 9) , (5) 
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Figure 2: Phase diagrams of the model (1) for d = q = 3 (a), and for d = 2, q = 5 (b). The 
correlation length is infinite along PC2 and at C4; q phases coexist at generic points along the 
dashed lines, q + 1 at the points C3 and C5. 



and measure the deviation from critical temperature and the magnetic field, respectively. 

For 2 < q < q c the transition along the line joining C q to C q ~\ in Fig. 1 is continuous and 
the scaling action (4) with r = 0, h < describes a massless flow from Aq FT to A.q F ^ . The 
scaling limit around the infrared fixed point C q -\ is described by the action 



Am = A 



(3-1) 

CFT 



d d xe{x) + A 



(9-1) 



d d Xi 



where f ~ m d ~ x ^ ' is proportional to T — T c and A ~ m d ~ X <t' is proportional to 1/H. All 
the operators in the r.h.s. of (6) are S'g-i-invariant: e is relevant, while <ft is the most relevant 
of the infinitely many irrelevant operators (dots) which specify the massless flow at r = 0. 

In d = 3 the condition 2 < q < q c is not satisfied for integer values of q and the transition 
at H < is first order for q > 3. For (/ = 3 the correlation length is finite at C q and infinite 
at C q -\. Since the latter is an Ising critical point for which e is the only symmetry-preserving 
relevant operator, C2 must be an infrared fixed point whose scaling region is described by the 
action (6). The ultraviolet endpoint of the massless flow at f = must be a fixed point P 
located on the transition line (Fig. 2a). The nature of the transition then requires that P is an 
Ising tricritical point 4 . 

In d = 2 the field theoretical description can rely on a number of exact results. Baxter solved 
the zero- field g-state Potts model at T = T c on the square lattice and found that q c = 4 [7] . The 
scaling limit of the critical point up to q c was later identified [13] to correspond to the conformal 
field theory with central charge [14] 
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(x ) + ... 



(9-1) 



(0) 



1 



t(t + l) 

where the parameter t is related to q by the formula 

TT(t-l) 



y/g = 2 sin 



2(t + l) 



(7) 



(8) 



4 Since Ising tricriticality is described by a $ 6 Landau- Ginzburg potential for which d = 3 is the upper critical 
dimension, P is a Gaussian fixed point. 
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The scaling dimensions of the leading thermal and magnetic operators coincide with those of 
the operators 02,1 an d 4>{t-i)/2,(t+i)/2 m the conformal theory, and read [13, 15] 

2 V t J 8t(t + 1) V ; 

The action (4) with d = 2 describes the scaling region around C q for q < 4 and is known to 
correspond to an integrable quantum field theory for h = [9] (see next section). For r = 0, 
h < and 2 < q < 4 it describes a massless flow between an ultraviolet fixed point with central 
charge c(q) and an infrared fixed point with central charge c(q — 1). Around this latter fixed 
point we can use the action (6) with the irrelevant operator eft identified with the operator 
of the conformal classification. Its scaling dimension at an Sq-invariant critical point is 

r{q) - ( - 2 



Xf=2\l + -j . (10) 

A field theory description around C g _i is still possible at q = 5. Since X^' = 2, (f> is in this case 
a marginal operator and the scaling region around C4 is described by the action 

A ( c ) FT -fjd 2 xe{x) + \jd 2 x(t){x). (11) 

Since the two-dimensional 4-state Potts model does not admit a tricritical point 5 [16, 17], (fi acts 
as a marginally relevant perturbation and the transition is first order (Fig. 2b). The action (11) 
is integrable also for t = [18, 19] and can be used to describe exactly the transition close to 
C 4 . 



3 Exact scaling solution in zero field 

In the two-dimensional case, to which we restrict our attention from now on, the action (4) is 
integrable for h = and the solution can be found in the form of an exact, elastic and factorized 
S"-matrix for the relativistic particles of the associated (1 + l)-dimensional theory [20]. The 
scattering theories above and below T c must describe two different physical situations and, at 
the same time, must reflect the existence of a duality transformation [5, 6] relating the ordered 
and disordered phases. 



3.1 Ordered phase 

The S-matrix in the case of spontaneously broken symmetry (r < 0) was determined 6 in [9]. 
The q ferromagnetic ground states correspond in the field theory to degenerate vacua labelled 

5 A phase diagram analogous to that of Fig. 2a, with P corresponding to a (q— l)-Potts tricritical point, should 
be obtained for 4 < q < 5. P tends to C q as q — * 4 and to C q -\ as q — > 5. 

6 The first exact S-matrix for ^2,1 -perturbed conformal field theories was determined in [21]. The relation 
between this solution, which does not exploit S q invariance and relies on a different particle basis, and the one of 
[9] is explained in [22]. 
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by an index a = 1,2, ... ,q. The elementary excitations are then provided by kinks 7 K a p(9) 
interpolating between the vacua a and (3 (a / /?). The space of asymptotic states consists of 
multi-kink configurations of the type K aoai (0i)K aiOl2 (62) ■ ■ ■ K an ian (9 n ) (aj / aj+i) interpo- 
lating between the vacua olq and a n . As a consequence of S^-invariance, all the n-kink states 
fall into two topological sectors: the neutral sector, corresponding to Qo = a n , and the charged 
sector, corresponding to ao / a„. 

Integrability implies that the scattering processes are completely elastic and factorised into 
the product of two-kink interactions. An outgoing two-kink state can only differ from the ingoing 
one by the vacuum state between the kinks. Hence, the two-kink scattering can formally be 
described through the Faddeev-Zamolodchikov commutation relation 

K^ejKtffa) = £ ^(^12)^(^)^/3(^1), (12) 

where #12 = 9i — 9 2 , and S^{9\ 2 ) denote the two-body scattering amplitudes (Fig. 3a). S q 
invariance reduces to four the number of independent amplitudes, two for the charged and two 
for the neutral topological sector 

K^ie^K^) = So(evi)Y, K a8(P2)K6l3(8l) + ft ^12)^7(^)^7/3(^1) , «^ 

5^7 

Ka^K^) = S 2 (9 12 ) ^ K a s(9 2 )K 8 a(0i) + S 3 {e 12 )K al {9 2 )K^ a {9x) . (13) 

5^7 

Using the commutation relation (12) twice one obtains the unitarity constraint 

E S Z(8)S £ a U-°) = > (14) 
which amounts to the set of equations 

(q - 3)S o (0)£o(-0) + ft(0)ft(-0) = 1 , (15) 

(g - 4)S (9)S (-9) + S (9)S 1 (-9) + ft(0)S o (-0) = , (16) 

(q - 2)S 2 {9)S 2 (-9) + S 3 (9)S 3 (-9) = 1 , (17) 

(q - 3)S 2 (9)S 2 (-9) + S 3 (9)S 2 (-9) + S 2 (9)S 3 (-9) = . (18) 

Crossing symmetry provides the relations 

S (9) = S (iir - 9) (19) 
Si(0) = S 2 (iir-9) (20) 
S 3 (9) = S 3 (iir-9) . (21) 

Using these constraints together with the Yang-Baxter and bootstrap equations (that we do not 
reproduce here) the following expressions for the four elementary amplitudes were determined 



7 We parameterise on-shell momenta as p M = (m cosh 8, m sinh#), m being the mass of the kink. 
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in Ref. [9] 



So(0) = ^^(g-™) . n(™) t (22) 

y ' sinhA(0- sinhA(#- f ) \i* J 



sin ^ sinh X(8-in) f\9\ 
sin^ smh\(e-^f) J 



sin 2 ^ sinhA6» f X6\ 
UG)= . S . 3 , w „ - , (24) 



sin^ sinh A (0- f ) \ 
A0\ 

in / 



7 7T 



sin A7T „ f X9\ 

^) = — n( — ) , (25) 



where A is related to g as 



and 



V^ = 2sin^, (26) 



\in)~ smhX(9-in) 6 ' 10 

^ f°°dx sinh i) -sinh f) x6 

AO) = / 2 V . 1 a; , r-f^ ^ sinh — . 28 

Jo x sinh cosh | m 

The above solution is well defined for real values of A, and one sees that (26) implies < q < 4. 
The pure scaling Potts model in this range of q is identified with the values of A in the range 
gome from 8 to 3/2. 

The function Il(A0/i7r) is free of poles in the physical strip Im 9 G (0,7r) for q < 3 (i.e. A < 1). 
Hence, in this range of q the only poles of the scattering amplitudes in the physical strip are 
those located at 9 = 2iir/3 and 9 = in/ 3 and correspond to the appearance of the elementary 
kink itself as a bound state in the direct and crossed channel, respectively. 

For q > 3 (A > 1) a direct channel (positive residue) pole located at 8 = 2in, 

«-j(i4). (») 

enters the physical strip in the amplitudes 82(8) and 83(8). Such a pole must be accordingly 
associated to a (topologically neutral) kink-antikink bound state B with mass 

riis = 2m cos k . (30) 

The amplitudes Si(8) and S3 (8) exhibit the corresponding crossed channel (negative residue) 
pole at 8 = in — 2in. The amplitudes Skb(8) and Sbb(8) describing the kink-bound state 
scattering and the bound state self-interaction are determined by the bootstrap equations 

Skb(9) = (q~ 2)S 2 (0 - iK)S!(8 + in) + S 3 (8 - i K )S 3 (8 + in) , 

S BB (0) = S BK (8-iK)S B K(O + iK), (31) 



3 The range 3/2 < A < 3 corresponds to the thermal perturbation of the tricritical q-state Potts model [9]. 
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(a) (b) 

Figure 3: Pictorial representation of the kink scattering amplitudes (a), and of the high- 
temperature amplitudes S^~^ k (b). 

and read 

Sbk(0) = ti-K/ir{Q)h/3-K/ir{Q) > ( 32 ) 

Sbb(0) = t 2/3 (9)t^ 2K/7T (e)t 2/ ^ 2K/7r (e) , (33) 

in terms of the functions 

, N tanh h(6 + ina) , ( 

t a (9) = Q '- . (34) 

tanh 2 (0 — iira) 

The poles located at = i(ir — k) in Sbk an d at 9 = 2ir/3 in Sbb are bound state poles 
corresponding to K and B, respectively. 

It has been shown [23] that the remaining poles in the amplitudes Skb and Sbb are associ- 
ated to multi-scattering processes rather than to new particles. Hence, the elementary kinks and 
their neutral bound state B are the only particles in the spectrum of the field theory describing 
the scaling zero-field Potts model below the critical temperature. 

3.2 Disordered phase 

We now show how the amplitudes (22)-(25) determine also the scattering theory in the disor- 
dered phase. 

Above the critical temperature there is a unique ground state and the excitations of the 
scaling limit must be ordinary particles rather than kinks. Observing that the group Z q of 
cyclic permutations is a subgroup of S q , we can take these particles to have a well defined Z q 
charge. The simplest possibility is then to conjecture a basis of elementary excitations Ai t (9), 
k = 1, . . . , q — 1, each carrying k units of Z q charge. A multi-particle state A^iOi) . . . Ak n {9 n ) 
will carry a charge k\ + . . . + k n (modg). The integrable scattering theory is characterized by 
the two-particle amplitudes defined through the Faddeev-Zamolodchikov algebra (Fig. 3b) 

A^Ap-M = £ Sljli (e 12 )A p _ j (9 2 )A j (9 1 ) , (35) 

where all the particle indices are taken mod q. The neutral channel corresponds to p = q, while 
p = 1, . . . q — 1 yields charged channels. Full S q invariance is then recovered requiring that the 
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interaction does not distinguish between the charged channels; within the neutral channel and 
the charged channel, the need to have well defined crossing properties forces to distiguish the 
case j = k from the case j 7^ k. This leaves us with a total number of four different amplitudes 
corresponding to 

4k(0i)iVfc(fc) = £o(M Yl A j(.h) A P-j(9i) + Si(M,4 fc (0 2 )A,_ fc (0i) > v^q 

A k {Q x )A q _ h {Q 2 ) = S 2 (9 l2 ) + S 3 (0i 2 )A k {0 2 ) A q _ k {0 x ) . (36) 

These identifications require that all the particles A k have the same mass m. 
The unitarity equations 

¥p 

are obtained iterating (35) and take a form identical to (15)-(18) with the substitution Si — > <Sj. 
The crossing relations coincide with (19)-(21) under the same substitution. 

This correspondence between the scattering theories above and below the critical temperature 
extends to the factorization and bootstrap equations and simply expresses the fact that the two 
phases are related by duality. The elementary excitations of the disordered phase are particles 
A k , k = 1, . . . , q — 1, with the same mass of the kinks of the low-temperature phase and whose 
scattering is expressed in terms of the same amplitudes which specify the kink S-matrix, i.e. 

m = m, Si(6) = Si(6), i = 0,1,2,3. (38) 

For 3 < q < 4, the high-temperature theory contains a neutral bound state B with the mass 
(30) and whose scattering is specified by amplitudes Sba u and Sbb coinciding with (32) and 
(33). 

Duality allows to compute correlation functions above and below T c within the form factor 
approach which relies on the knowledge of the S-matrix. This was done in [11] using the kink 
S"-matrix. It can be checked that the same results are obtained through the high-temperature 
scattering theory. 



4 Particle spectrum in weak magnetic field 

The action (4) describes the renormalization group trajectories flowing out of the fixed point 
located at the origin of the r-h plane. Such trajectories can be labelled by the dimensionless 
parameters 

^± = (± h ){d-X e )/(d-X a ) ' ( 39 ) 

where the upper and lower signs are used for h > and h < 0, respectively, in such a way that 
•q + parameterizes the trajectories in the upper half-plane and ?]_ those in the lower half-plane; 
the two trajectories at h = correspond to rj± = +00 and 77-1- = —00 (see Fig. 4). 
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Figure 4: Some renormalization group trajectories associated to the action (4). 

In this section we discuss the evolution of the mass spectrum in d = 2 starting from the case 
of small magnetic field for which perturbation theory around integrable quantum field theories 
[10] can be used. 

4.1 Weak field above critical temperature 

For small h at r > (i.e. for rj± — > +oo) the corrections to the mass spectrum are determined by 
the matrix elements of the magnetic operator cr q (x) on the asymptotic states of the unperturbed 
(rj± = +oo) theory. At leading order, the correction to the mass matrix is [10] (Fig. 5a) 

(5m\ k ~ -2^(0)1^(0)1^(0)), j,k = l,...,q-l (40) 

(here and below the matrix element on the particles A k are intended at h = 0). 
Let us introduce the operators 

T k (x) = Y,u- ak a a (x), k = l,...,q-l (41) 

a=l 

where uj n = exp(2i7r/n). The generator Q q of cyclic permutations acts as 

n q a a (x) = a a+l (mod q) (x) , n q r k (x) = uj* T k (x) , (42) 

showing that r k carries k units of Z q charge and can be taken as interpolating operator of particle 
A k at h = 0. Using 

n 

E^=«4,n (43) 

and (3), one also have 

1 9-1 

Vq(x) = - J2 T k(%) ■ (44) 
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Figure 5: First order mass corrections in weak magnetic field above critical temperature (a), 
and below (b). 

Denoting 

A$ k = (Aj(0)\*(0)\M0)), (45) 
and taking into account conservation of Z q charge as well as full S q symmetry at h = 0, we have 

Mj; k = f q 5 jMl{modq) (46) 

and 

M fi> = ^T, M j!k = %( 1 -Si*) ( 47 ) 
q i=i q 

for the matrix which determines the mass corrections (40). Diagonalization of this matrix gives 
for the spectrum at r)± — > +oo a particle 

1 q ^ 

Ao = £ A k ( 48 ) 



with square mass 



and a degenerate multiplet 



ml ~ m 2 - 2f q - — - h , (49) 



A' k = ^=(A k -A q - 1 ), k = l,...,q-2 (50) 



with square mass 



f 

m' 2 ~ m 2 + 2 J -± h . (51) 

Comparison with (44) shows that Aq is interpolated by a q and is a singlet of the S q -\ symmetry 
surviving at h ^ 0; suitable linear combinations of the A' k yield a multiplet in which each of the 
q — 2 particles carries a definite (non-zero) Z q _\ charge. 

For 3 < q < 4 the theory also contains the additional S'g-i-singlet B with a mass correction 
with respect to (30) given by 

Sm 2 B ~ -2h (B(0)|<7,(0)|B(0)) . (52) 

This matrix element, as well as f q above, can be obtained from the form factor results of [11] 
(for q = 4 see also [8]). Notice that the first order correction to the mass of Aq vanishes for 
q = 2; the leading correction in this case is of order h 2 and was computed in [24]. 
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4.2 Weak field below critical temperature 

A small magnetic field acting on the sites with color q affects the ground state degeneracy at 
r < 0. Denoting by \0 a ), a = 1, . . . , q the ferromagnetic vacua, S q symmetry gives at r]± = — oo 

v 

(cr 7 ) Q = (0 a \a y {x)\0 a ) = — — j- (q5 Jta - 1) , (53) 

with v positive. At first order in h, the energy density difference between a vacuum \0 a ^ q ) and 
|0 g ) is then 

AS = S£ a - d£ q ~ -h ({a q ) a - (a q ) q ) = -^L h , (54) 

so that |0 g ) is the unique true vacuum at h > 0, and the unique false vacuum at h < 0. 

Since no finite-energy topological excitation can begin or end on a false vacuum, the space 
of asymptotic states of the theory does not contain kinks for h > 0. For large and negative, 
instead, the kinks K a p with a, (3 7^ q survive as the elementary excitations of the theory. The 
first order correction to their mass is (Fig 5b) 

Sm 2 Kap ^ -2h <A>(0)|a 9 (0)|i^ Q (0)) corm = -2hF^J^) , (55) 

where 

Kpai ! ~ h) = (v«W q (Q)\K a p(ei)K Pa (e 2 )) , a + p (56) 

is the two-kink form factor at rj± = — oo [11]. In general this function has an annihilation pole 
whose residue 

-iRes e= in F a p a (0) = {° q )a - W q ) f3 = -—^ (S a>q - 6p >q ) (57) 

vanishes precisely when both a and f3 differ from q, giving a finite mass correction 9 (55) (which 
does not depend on a,j3 = 1, . . . , q — 1). The divergence of (55) when a or (3 equal q simply 
reflects the decoupling of a kink interpolating between vacua which are no longer degenerate. 

Although for h positive |0 g ) is the only true vacuum and no single-kink state survives, n- 
kink states beginning and ending on the true vacuum, i.e. K qai (6i)K aia2 (02) . ■ ■ K an _ iq (6 n ), do 
not decouple. Here we consider all the intermediate vacua to be false, since insertion of a true 
intermediate vacuum would simply amount to breaking the sequence into two states of the first 
type. Then we have n < q. 

In the non-relativistic limit valid for small rapidities, the total energy of such a configuration 
consist of the rest mass term nm, the kinetic term, a contribution coming from kink interaction 
which decays exponentially with interkink distance, and the false vacuum contribution 

n-l 

V(xi, ...,x n ) = A£ y~](x i+ i - Xi) = (x n - xi) AS , (58) 

i=l 

9 In this case (55) is finite irrespectively of the sign of h. We should recall, however, that for ft > we are 
computing a mass gap above a false vacuum whose energy separation from the true vacuum in a system of spatial 
size L is A£ L. Hence, all the single-kink states decouple in the infinite system with a positive, however small, 
magnetic field. 
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q a q 



(a) 



• • • 

q a 1 a 2 q 

(b) 



Figure 6: Spatial kink configurations which are confined into mesons (a) and baryons (b) by a 
positive magnetic field. 

where x\ < X2 < ■ ■ ■ < x n are the spatial positions of the kinks. The positive linear potential (58) 
confines the n kinks into a finite spatial interval and prevents the observation of isolated kinks. In 
this sense the kinks of the r/± = — oo theory play the role of "quarks" at h > 0. The asymptotic 
particles are instead the topologically neutral bound states produced by the confinement of the 
n-kink state. When rj + — > — oo the confining potential is extremely shallow and the kinks are 
very loosely bound; the average interkink distance is large and kink interaction is negligible in 
first approximation. The n-kink bound states form an infinite tower of levels which are dense 
above the value nm as rj + — > — oo. It is natural to call "mesons" the n = 2 (kink-antikink) 
bound states (Fig. 6a), and "baryons" the bound states with n = 3 (Fig. 6b). Recalling the 
conditions n < q and q < q c = 4, we have that the mesons can occur for q = 2, 3, 4, and the 
baryons for q = 3, 4; tetraquark confined states are allowed for q = 4. 

These particles organize themselves into multiplets of the residual S q -i symmetry. The 
number of different n-kink sequences coincides with the possible ways of coloring the intermediate 
vacua. The q — 1 mesonic sequences can be combined into the states 

q-l 

4 J) (0) ~ E K qa (9)K aq (-9) , k = 0, 1, . . . , q - 2 (59) 

a=l 

with j = 1,2,... labelling in order of increasing energy the levels originated by the confinement 
of the kink-antikink superposition. Since 

fVi^^-i^, (60) 

the mesons (59) are eigenstates under cyclic permutations of the first q — 1 colors with Z q ^\ 
charge k. In the non-relativistic limit valid for the lowest levels in weak field, we can think 
of the kink and antikink inside a meson as experiencing elastic reflection on the walls of the 
confining potential and elastic scattering among themselves. As the effect of the latter, the 
intermediate vacuum can either remain unchanged with a probability amplitude S3 (2$), or 
switch to a different color with a probability amplitude E2(26>) which, by 5,j_i-invariance, does 
not depend on the new color 10 . The superpositions in the right hand side of (59) are eigenstates 
of the S-matrix with scattering amplitudes 

£ 3 (2#) + [(q- l)S k , - 1]E 2 (20) . (61) 

10 At rj± — —00 the scattering amplitudes £2 and £3 coincide with (24) and (25), respectively. Corrections in 
weak field are determined by the matrix elements of the magnetic operator a q [10] . 
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Figure 7: Seven-kink matrix element responsible for the decay of baryons above threshold into 
two mesons in positive magnetic field. 



Hence, the quark interaction is different for the neutral mesons {k = 0) and for the charged 
mesons (k ^ 0). For fixed j this will lead to a siglet 7Tq and a multiplet tt^ , . . . ,^ql 2 with 
energies which differ very slightly in weak field. 

Only a finite number of these confined states are stable. For q = 2,3 all those states with 
energy larger than twice the mass of the lightest mesons (i.e. larger than a value close to 4m for 
weak field) lie in the continuum and can decay. This means that for r/ + sufficiently large and 
negative we certainly have stable mesons in the energy interval (2m, 4m) for q = 2, 3, and stable 
baryons in the energy interval (3m, 4m) for q = 3. At q = 4 the decay thresholds are lowered by 
the presence of the neutral particle B with unperturbed 11 mass mg = y/3m. Hence the stability 
intervals in weak field are (2m, 2\/3m) for the neutral mesons, (2m, (2 + y/3)m) for the charged 
mesons, (3m, 2\/3m) for the neutral baryons and (3m, (2 + \/3)m) for the charged baryons; as 
for the tetraquark states, they all lie in the continuum and are expected to be unstable. 

The number of 3-kink sequences giving rise to baryons is (q — l)(q — 2). For q = 3 we then 
have the two series of baryons 



P { £\e) ~ K 31 (e 1 )K 12 {e 2 )K 23 {9 3 ) ± K 32 (e 1 )K 21 (e 2 )K 13 (e 3 



(62) 



with even or odd parity with respect the residual Z 2 symmetry which interchanges the colors 
1 and 2. For q = 4 the residual S 3 symmetry can be seen as the product of the group Z 3 of 
cyclic permutations times the topological charge conjugation Ct which transforms the kink K a g, 
a, j3 7^ 4, into its antikink Kp a . Then the six series of baryonic states 



-ka 



a=l 



are ^-eigenstates: 



^4, a (6 l l)^ Qia+1(mo d3)( 6 '2)^ Q+1(m od3),4( 6, 3) ± 
K A,a+l (mod 3) (01 ) K a+l (mod 3),a ^) K »A (^) 



U 3Pk,± - W 3 Pk,± 



fc = 0,1,2 (63) 
(64) 



L The first order correction to m% coincides with (52). 
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For a given j they give rise to two singlets (k = 0) and two doublets {k ^ 0) with opposite 
C^-parity. 

As already said, a small negative magnetic field confines only the kinks K a p with a or (3 
equal q. Since now the false vacuum is unique, the only confined states are mesons of type 
K aq K qa . For q > 2, however, these can decay into two asymptotic kinks through expansion of 
bubbles of true vacuum into the false vacuum. 

5 Evolution of the spectrum with temperature and magnetic 
field 

5.1 Positive magnetic field 

We discussed in the previous section how the particle spectrum of the scaling two-dimensional 
q-state Potts model changes when a small magnetic field acting on the sites with color q is 
switched on. For positive field, below critical temperature we have a complete confinement of 
kinks and the generation of a dense spectrum of mesons and (for q > 2) baryons; the number of 
such particles which are stable tends to infinity as 77+ — > —00. Above critical temperature, on 
the other hand, the effect of the magnetic field is much less dramatic, simply amounting to a 
partial removal of the degeneracy of the mass spectrum of the h = particles. 

It must be possible to interpolate continuously between these two limiting cases following 
the evolution of the spectrum as rj + grows from —00 to +00 on the plane of Fig. 4. The simplest 
scenario is that, as r/+ increases, more and more mesons and baryons cross the decay thresholds 
and become unstable. This process of depletion of the spectrum of stable excitations would 
continue until the only stable particles surviving as r/ + — > +00 are those of the high-temperature 
theory in zero field. 

We expect that, during the evolution of the spectrum as a function of the parameters rj±, 
energy levels corresponding to states originating from the confinement of a same number of 
kinks do not cross for h 7^ 0. Indeed, the additional degeneracy at a crossing point of this 
nature should normally be related to a symmetry enhancement. In the field theory (4) the 
only symmetry enhancement (from S q -i to S q ) occurs at h = 0. If we add to this that the 
baryons should decay more easily than the lightest mesons, it seems natural to expect that 
the particles surviving in the limit r/ + — > +00 should be identified with the lightest among the 
mesons produced by kink confinement at very large and negative. In particular, the lightest 
meson multiplet , . . . , 7r^_ 2 should evolve for increasing i] + into the multiplet (50). As for 
the particle which evolves into the singlet (48), it should be identified with 7^ for q = 2, 3. 
At q = 4 the zero-field spectrum also includes the neutral particle B with mass v3m, which is 
the lightest particle for rj + — > — 00. If one assumes that it does not cross mesonic levels, then 
it should evolve into the neutral particle (48) for rj + — > +00; then the meson tt^ could evolve 
into the particle with mass \/3m in the same limit. 

According to this scenario, all the mesons (59) with j > 1 and all the baryons must have 
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become unstable by the time 77+ approaches +00. For each of this particles there should exist 
a finite "critical" value of t/ + for which they reach the lowest decay threshold compatible with 
their charge, and above which they become unstable. Such critical value is expected to decrease 
as the mass of the particle increases, so that the "critical trajectories" accumulate in the limit 
77+ — > —00. Notice that, due to the non- locality of the magnetic operator with respect to the 
kinks, the magnetic term of the action has infinitely many matrix elements on kinks which are 
non-zero in zero-field. In particular, the vertex responsible for the decay of the baryons which 
reach the two- meson threshold is shown in (Fig. 7). 

Also in view of considerations to be made below about the spectrum evolution at h < 0, we 
expect the quantity f q determining the first order mass corrections (49) and (51) to be positive 12 . 
Together with the previous identifications, this leads to the expectation that for any h > the 
different quark interactions inside the neutral and the charged mesons (see (61)) induce a mass 
splitting between n^, . . . , 7T^ 2 an d which is positive for q = 3. If the previous speculations 
about g = 4 should turn out to be correct, the mass splitting would be negative in this case. 

When specialized to q = 2 this scenario coincides with that originally proposed for the 
scaling Ising model by McCoy and Wu [3]. In this case the pattern is simplified by the absence 
of baryons and charged mesons, as well as by the absence of interaction at h = 0. Moreover, 
integrability at r = [25] also allows analytic investigation for strong magnetic field. Several 
studies, both analytic and numerical, have confirmed the McCoy- Wu scenario and provide us 
with a detailed description of the mass spectrum of Ising field theory in the full r-h plane 
[3, 25, 10, 26, 24, 27, 28, 29, 30, 31] (see also [32] for an introductory review of Ising field 
theory). 

5.2 Negative magnetic field 

We now discuss the evolution of the spectrum as a function of 77 _ for q = 3, 4 (for q = 2 the 
spectrum does not depend on the sign of the magnetic field). We saw in the previous section 
that the elementary excitations for t?_ — > — 00 are the kinks K a p interpolating between the 
q — 1 residual vacua. For t?_ — > +00, on the other hand, due to the negative sign of h in (49) 
and (51), the lightest excitations are the charged particles that in the previous subsection we 
identified with the mesons , . . . , tt^ 2 - Passing from these topologically neutral particles to 
the topologically charged kinks requires a phase transition at an intermediate value of rj- . We 
already know that, for q = 3,4, a second order phase transition takes place at r/_ = 13 . The 
mass of the kinks decreases as r/_ goes from —00 to zero, and vanishes at r\- = 0. Similarly, 
the mass of the mesons tt^ , . . . , 7r^ 2 decreases as r/_ goes from +00 to zero, and vanishes at 
7/_ = 0. Deconfined kinks exist only for r/_ 6 [— 00, 0); for r/_ £ (0, +00] and for all positive 
values of h the spectrum can be seen as consisting of kink bound states. The trajectory t?_ = is 
a massless flow between the S^-invariant ultraviolet fixed point at h = and the S g _i-invariant 

12 This quantity can be exactly computed within the integrable field theory at h — (see [11]). 
13 An analogous transition from ordinary particles to kinks, involving the spontaneous breaking of the symmetry 
S q rather than 5*5-1, takes place when going through r = at h — 0. 
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Figure 8: Conjectured qualitative evolution of the mass spectrum for q = 3. The left-hand 
half corresponds to the evolution in rj + , the right-hand half to the evolution in ?]_ . The values 
r]± = +00, as well as the values ij± = —00, describe the same renormalization group trajectory. 
Only three lightest mesons and the lightest neutral baryon are shown for positive magnetic field. 
The dashed line is the lowest decay threshold (twice the mass of the lightest particle). Particles 
whose mass reaches the threshold become unstable. Deconfined kinks exist for negative ry_ only. 



infrared fixed point at h = —00. 

The ratio between the mass of the neutral meson ir^ and that of the charged mesons 
tt^\ . . . , 7Tg_ 2 , which is 1 at r?_ = +00, increases as r\- decreases until it reaches the value 2 for 
some positive value of 7]__. Beyond this point the opening of the decay channel ir^ — > n^ni 1 ^ 
makes the neutral meson unstable. The second neutral particle present at q = 4 becomes 
unstable in the same way. We then expect that the charged mesons are the only stable particles 
surviving for r}- sufficiently small and positive. Similarly for the kinks at sufficiently small 
and negative. 

The conjectured evolution of the mass spectrum with the parameters 77+ and ?y_ is illustrated 
in Fig. 8 for q = 3. 



6 Conclusion 

In this paper we considered the field theory describing the scaling limit of the two-dimensional 
g-state Potts model in a magnetic field acting on one of the q colors. This field breaks the 
S q symmetry of color permutations down to S q ~i and, for 2 < q < q c = 4 allows for an 
extended region in the plane of temperature and magnetic field in which the quarks (kinks 
interpolating between degenerate vacua) are deconfined. If the analysis is extended to a more 
general magnetic term Yla^ 1 h a &a, a renormalization group trajectory corresponding to generic 
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values of the components h a will possess no internal symmetry and will contain only particles 
made of confined quarks with unequal masses. Of course, suitable relations among the magnetic 
parameters identify regions in parameter space with a residual S q -i or (for q = 4) Sq-% symmetry, 
which in turn contain phases with deconfined quarks. 

We saw how form factor perturbation theory allows to compute mass corrections in weak 
field above critical temperature, and below critical temperature for the quarks in the deconfined 
phase. As for the mass spectrum of the particles originating from confinement, the Bethe- 
Salpeter approach proved to be remarkably effective for the Ising mesons [26, 29], but needs to 
be generalized to the case of quarks which interact already in absence of field in order to deal 
with q ^ 2. Hopefully, it will become possible to study also the baryonic spectrum along similar 
lines. 

Numerical methods will be eventually needed for a quantitative study of the mass spectrum in 
strong field. Particularly promising in this respect seems the truncated conformal space approach 
[33] , which amounts to the numerical diagonalization of the Hamiltonian on a finite-dimensional 
subspace of the conformal basis of states of the ultraviolet fixed point. This approach has been 
successfully used for the Ising case 14 [10, 26, 29, 34] and can be used also for q ^ 2. 

Acknowledgments. The work of G.D. is partially supported by the ESF grant INSTANS 
and by the MUR project "Quantum field theory and statistical mechanics in low dimensions". 



14 In [26, 29] the free nature of the zero-field Ising model is exploited to diagonalize the Hamiltonian on a 
truncated basis of massive fermionic states. 
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